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1. Introduction 

Supersymmetric field theories exhibit many remarkable properties. Chief among these are the 
cancellations which occur between boson and fermion contributions in the perturbative calculation 
of physical quantities. These cancellations eliminate many of the divergences typical of quantum 
field theory and are at the heart of its use to solve the gauge hierarchy problem [[I]]. Additionally, 
supersymmetric versions of Yang Mills theory, while more tractable analytically than their non- 
supersymmetric counterparts, exhibit many of the same features such as confinement and chiral 
symmetry breaking [0]. Super Yang-Mills theories with extended supersymmetry are also conjec- 
tured to be dual to various supergravity theories in the limit of a large number of colors and have 
been proposed as a regularization for M-theory [||. [|] . 

Since low energy physics is manifestly not supersymmetric it is necessary that this symmetry 
be broken at some energy scale. A set of non-renormalization theorems ensure that if SUSY is not 
broken at tree level then it cannot be broken in any finite order of perturbation theory see eg. [fi]]. 
Thus the general expectation is that any breaking should occur non-perturbatively. The lattice 
furnishes the only tool for a systematic investigation of non-perturbative effects in field theories 
and so significant effort has gone into formulating SUSY theories on the lattice - see ^ and the 
recent reviews 

Unfortunately, there are several barriers to such lattice formulations. Firstly, supersymmetry is 
a spacetime symmetry which is generically broken by the discretization procedure. In this it resem- 
bles Poincare invariance which is also not preserved in a lattice theory. However, unlike Poincare 
invariance there is usually no SUSY analog of the discrete translation and cubic rotation groups 
which are left unbroken on the lattice. In the latter case the existence of these remaining discrete 
symmetries is sufficient to prohibit the appearance of relevant operators in the long wavelength lat- 
tice effective action which violate the full symmetry group. This ensures that Poincare invariance 
is achieved automatically without fine tuning in the continuum limit. Since generic latticizations of 
supersymmetric theories do not have this property their effective actions typically contain relevant 
supersymmetry breaking interactions. To achieve a supersymmetric continuum limit then requires 
fine tuning the bare lattice couplings of all these SUSY violating terms - typically a very difficult 
proposition. 

Secondly, supersymmetric theories necessarily involve fermionic fields which generically suf- 
fer from doubling problems when we attempt to define them on the lattice. The presence of extra 
fermionic modes furnishes yet another source of supersymmetry breaking since typically they are 
not paired with corresponding bosonic states. Furthermore, most methods of eliminating the extra 
fermionic modes serve to break supersymmetry also. 

Recently there has been renewed interest in this problem stemming from the realization that in 
certain classes of theory it may be possible to preserve at least some of the supersymmetry exactly 
at finite lattice spacing. It is hoped that this residual supersymmetry may protect the lattice theory 
from the dangerous SUSY-violating radiative corrections that leading to fine tuning. In the case of 
extended supersymmetry two approaches have been followed 1 ; in the first, pioneered by Kaplan 
and collaborators, the lattice theory is constructed by orbifolding a certain supersymmetric matrix 
model and then extracting the lattice theory by expansion around some vacuum state [[K], 11, 12, 



'For examples of ,jV = 1 models with exact SUSY see the recent work |^] and |^] 
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13]. The second approach relies on reformulating the supersymmetric theory in terms of a new set 



of variables - the twisted fields. In this procedure a scalar nilpotent supercharge is exposed and 
it is the algebra of this charge that one may hope to preserve under discretization [14], [15], [|k]]. 
This approach was initially used to construct and study lattice formulations of a variety of low 
dimensional theories without gauge symmetry [17, |TJ|, [l^, 2C, 21 j. Subsequently it was extended 
to case of Yang-Mills theories by Sugino [22, 23] although that construction was confined to low 
dimensional theories. In this review we will confine ourselves to a discussion of a new geometrical 
discretization of the twisted theory which accommodates both two and four dimensional super 
Yang Mills theories. 

The basic idea in this approach is decompose the fermion fields as a set of antisymmetric 
tensor fields which are subsequently interpreted as components of single Dirac-Kahler field. This 
same procedure when applied to the supercharges exposes a scalar nilpotent supercharge Q. Fur- 
thermore, the action of the supersymmetric field theory can then be written as the <2-variation of 
some function. Provided we can then maintain the nilpotent property of Q under discretization it 
is straightforward to write down a lattice action in terms of these twisted fields which is exactly 
invariant under Q at finite lattice spacing. As a bonus we will find that that this twisted formulation 
of the theory, being written entirely in geometrical terms, can be discretized without encountering 
spectrum doubling. 

We first describe the nature of the twist in two dimensions and explain how it is related to 
the use of Dirac-Kahler fermions. We then give a detailed discussion of continuum jV = 2 super 
Yang-Mills in two dimensions which furnishes a simple example of the twisted approach. Our 
discretization prescription is then introduced and a lattice action for this theory constructed. The 
latter is local, free of spectrum doubling, gauge invariant and preserves a single (twisted) super- 
symmetry at finite lattice spacing. It is naturally formulated in terms of complex fields and we 
discuss how this lattice model may be used to target the true continuum theory. Preliminary results 
are then shown coming from a RHMC simulation of this lattice model which indeed show an exact 
<2-symmetry emerging at large coupling j8 . The technique can be generalized to provide a lattice 
action for jV = 4 super Yang-Mills in four dimensions and this is then discussed in some detail. 



2. The 2D Twist 

Consider a theory with jV = 2 supersymmetry in two dimensional Euclidean space. Such 
a theory contains two Majorana supercharges q ! a transforming under the global symmetry group 
SO (2) x SOr(2) where the subscript corresponds to two dimensional rotations and the superscript 
describes the behavior under the R-symmetry corresponding to rotating the two Majorana fields 
into one another. 



The basic idea of twisting which goes back to Witten [ 24 ] is to introduce a new rotation group 



SO(2)' = diagonal subgroup (50(2) x 50(2)*) (2.1) 

and to decompose all fields now as representations of this new rotation group. In simple terms 
what this means is that whenever I do a rotation in the base space by some angle I must do an 
equal rotation in the R-symmetry space. Thus it is equivalent to treating the two indices / and 
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a as equivalent. Thus the supercharges are to be interpreted as matrices in this twisted picture 

<fa <7a/3 (2-2) 
It is now natural to expand such a matrix on a basis of products of 2D gamma matrices 

q = QI + QnYn + QnYiYi (2.3) 

The fields (Q,Q^,Qn) w& called the twisted supercharges. The original SUSY algebra then im- 
plies a corresponding twisted algebra which takes the form 

{<?,<?}«£ =4^" (2.4) 

In components this reads 

{Q,Q} = {Qn,Qn} = {Q,Qn} = {Qn,Qv} = (2.5) 

{Q,Q^i}=P ll (2.6) 

{Qn,Q»} = -£nvPv (2-7) 

Notice that the scalar component of the supercharge matrix is nilpotent as previously advertised. 
It is also important to realize that the twisted superalgebra implies that the momentum P is now 
the <2-variation of something - it is said to be Q-exact. This fact renders it plausible that the entire 
energy momentum tensor may be 2-exact in such twisted theories. Since 7^ v = this would 
imply that the entire action of the theory could be written in a <2-exact form S = QA [|27|]. This 
turns out to be true in essentially all the theories that have been studied. Specifically we will show 
that the action of ,jV = 2 super Yang-Mills theory in two dimensions can be written in this way 



[28]. Notice also that to match the four real supercharges of the original supersymmetric theory it 
is natural to take the twisted supercharges to be real also. This in turn implies that the supercharge 
matrix obeys a reality condition 

q = q* =q T m2D (2.8) 

If the presence of a nilpotent scalar fermionic charge and a corresponding Q-exact action all 
sound reminiscent of BRST gauge fixing you would be right. These twisted theories can be derived 



by gauge fixing an underlying topological symmetry [|29|]. If the gauge fixing function is chosen 
carefully the resultant theory can be untwisted to reveal an associated supersymmetric field theory. 
However, it must be stressed that this correspondence is not one to one. In the topological field 
theory the set of physical states must be restricted in the usual way to those annihilated by the 
scalar charge. In the language of the supersymmetric field theory such a constraint would project 
one to the vacuum state. This explains why such a theory contains no physically propagating 
degrees of freedom - it is a theory defined only on the moduli space of classical solutions. Our 
proposal is not to impose this physical state condition but merely to treat the twisting procedure as 
an exotic change of variables in the supersymmetric field theory. Certainly in flat space the physical 
content of the twisted and untwisted theories is the same. 

Let us turn now to a discussion of the fields in such a twisted theory. It should be clear that if 
the supercharges can be written as matrix so can the fermions 

x V=^I + Wn 7n +X12Y1Y2 (2-9) 
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Thus the twisted theory will not contain spinors but antisymmetric tensor fields. Indeed it is pos- 
sible to abstract these p-form fields and consider them as components of a geometrical object - a 



Dirac-Kahler field as first emphasized by Kawamoto and collaborators [ ]16| ] . 



¥=(|,lMi 2 ) (2-10) 

It is a remarkable fact that the original Dirac equation for 2 (Majorana) fermions is then equivalent 
to the geometrical Dirac-Kahler equation 



(d-d J( ) x ¥ = (2.11) 

Here d and d^ are the usual exterior derivative and its adjoint. Their action on a generic Dirac- 
Kahler field co is given by 

d(O = {0,d^f,d ll f v -d v f tl ,...) (2.12) 

and 

-tfa> = (f ,/;,..., 0). v (2.13) 
The Dirac-Kahler equation can be derived from an action 

S F = ( x ¥ f .(d-d^) x ¥) (2.14) 

where the dot product of two Dirac-Kahler fields A and B is defined by 

<A\B>= ftPx^g^A* '^Mp ( 2 - 15 ) 

Writing out these components explicitly and restricting ourselves to flat space we find 

Sf = ^vJ^J] + -xlydfr Yv] (2-16) 

This equivalence of the Dirac equation to the Dirac-Kahler equation holds equally when ordinary 
derivatives are replaced by gauge covariant derivatives. 

The usual reason why the Dirac equation is preferred over the Dirac-Kahler equation is that 
the latter naturally describes the propagation of 2 degenerate fermions in two dimensions. In a 
theory like QCD where the bare quark masses are not usually thought of as equal this appears 
problematic. However, a theory with jV = 2 supersymmetry in D = 2 dimensions naturally requires 
such degenerate fermionic states. Thus we have argued that for these special theories the twisting 
procedure we have described is actually exactly equivalent to the decomposition of the fermion 
fields in terms a single Dirac-Kahler field. Furthermore the equivalence between the Dirac equation 
and the Dirac-Kahler equation holds also in four dimensions where the latter equations describes 
the propagation of four degenerate (Majorana) fermions. Again, this is precisely the field content 
of jV = 4 super Yang-Mills theory and we should not be surprised then to find out that the latter 
theory when suitably twisted can be expressed in terms of the components of a single real Dirac- 
Kahler field [U, 



Finally any theory invariant under the <2-symmetry must contain superpartners for these fermionic 
fields. In two dimensions we then expect to find the following bosonic (commuting) fields 

@,Ap,B n ) (2.17) 
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The nilpotency of Q then requires that Q 2 = up to gauge transformations when acting on any of 
these component fields. If we interpret Bu as a multiplier field we see that in this way we have 
generated the bosonic field content of jY = 2 SYM in 2D if we interpret each field as taking values 
in the adjoint of some gauge group. We turn now to an explicit construction of the twisted form of 
this theory. 



3. Twisted JV = 2 SYM in D = 2 

In this section we render all these ideas concrete by writing down the twisted formulation of 
JV = 2 super Yang-Mills theory in two (continuum) dimensions. The construction will expose the 
scalar supersymmetric invariance explicitly and the Q-exact form of the action. We will also show 
that the resulting on-shell action can be rewritten in terms of the usual spinor degrees of freedom 
exhibiting the fact that (in flat space) the twisted theory is nothing but a change of variables in the 



original spinor based super Yang-Mills model. The twisted gauge fermion is written [|28|],[|22|] 

S = pQTr J d 2 x(^[<l>,$}+2xi2Fi2+ (3.1) 

+ XnBn + Y^D^) (3.2) 

where each field is in the adjoint of some gauge U(N) group X = £^ =1 X f T" and we employ 
antihermitian generators T". The action of Q on the fields is given by 





= Vn 


QWn 


= -Drf 




= 


Qln 


= #12 


QBn 




& 


= n 


On 


= [0,0] 



(3.3) 

Notice that Q 2 generates an (infinitessimal) gauge transformation parametrized by the field (p. Car- 
rying out the <2-variation and subsequently integrating out the multiplier field B\2 we find the 
on-shell action 

s = pTrJ d 2 x(hpM 2 -\n^>,n\-Fh (3-4) 

- D^D^-Xn[$,Xn] (3.5) 

- 2xi 2 {DiW2-D 2 Wi)-2Y fl D fl T]/2 (3.6) 

The gauge and scalar part of this action are easily recognized to be the usual bosonic sector of 
jY = 2 super Yang-Mills. If I choose the contour where 0° = (0 a )* it is clear that the bosonic 
action is real and positive semi-definite. The fermionic sector looks a little unconventional but by 
looking back at the previous discussion it can be seen that it merely corresponds to a Dirac-Kahler 
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decomposition of the usual spinor action corresponding to two Majorana fermions. To convince 
ourselves of the equivalence of this twisted formulation to the usual one it is convenient to pick a 
Euclidean chiral basis for the two dimensional /-matrices 




W=l 10 )' 
and to subsequently define a Dirac spinor A as 

(3.9) 



Vi - m 



It is then a straightforward exercise to show that the Dirac action X^y.DX reduces to the kinetic 
terms in the twisted action we have previously uncovered. Indeed the entire twisted action including 
the Yukawas can be rewritten in terms of this Dirac spinor and can be recognized as the usual 
fermionic sector of two dimensional JV = 2 super Yang-Mills. 

5 = A t Af(0)A (3.10) 

where 

M = r .D+<I±%,H-^[*,H (3.1D 

4. Discretization 

We will adopt this geometrical formulation as a useful starting off point for constructing a lat- 
tice theory. It is convenient from the point of view of supersymmetry as the scalar component of the 
twisted supersymmetry involves no derivatives and offers the hope, therefore, of being compatible 
with a lattice structure. Secondly, the use of Dirac-Kahler fermions is attractive as it is known how 
to discretize the Dirac-Kahler action without inducing spectrum doubling. Indeed, Dirac-Kahler 
fermions at the level of free field theory are known to be equivalent to staggered fermions. The 
usual degeneracy of staggered fermions now becomes a plus - it allows us to replicated the degen- 
erate flavors required by extended supersymmetry. The use of Dirac-Kahler fermions to formulate 



(Hamiltonian) lattice supersymmetric models was first proposed in [33]. 

Clearly, continuum scalar fields should map to fields defined on sites in a (hypercubic) lattice, 
continuum vector fields to fields defined on links and rank 2 antisymmetric fields to lattice fields 
living on plaquettes. Notice that all lattice p-cubes whose dimension is greater than zero come 
in two possible orientations corresponding to whether the defining vertices are written down as 
an even or odd permutation of some standard ordering. In the case of links this is just the usual 
statement that the link has two possible directions. It is then natural (and as we will see later 
essential) that this 2-fold orientation of the underlying p-cube can be associated with the existence 
of two independent lattice fields for each continuum field. These two fields can be represented 
by allowing the lattice fields X" to be complex numbers thus promoting the group from U (N) to 
GL(N,C). The lattice fields are also assigned gauge transformation properties which generalize the 
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usual transformation properties of site fields and links in lattice gauge theory [ 34 ] : 

f(x) -»• G(x)f ( x)G-\x) (4.1) 
f\i to G(x)/ P (*)<T *(* + M) (4-2) 
/ 12 (x) - G(i)/ 12 (i)G-'(i+1+2) (4.3) 

where G = is a lattice gauge transformation. Notice in the naive continuum limit where G ~ 
1 + these reduce to the usual commutator structure typical of adjoint fields. In general this leads 
to a point split or shifted commutator. For link fields this looks like 

[<K*)/ m (*)-//*M<K*+M)] (4-4) 

Notice that this commutator has a well-defined lattice gauge transformation property. Finally, we 
must replace the continuum gauge connection with a Wilson gauge link U^(x) = er*. Notice 
now that will be complex which implies that the gauge links are not unitary at this stage in the 
construction. 

Of course we have not yet finished. We need also to give a prescription for replacing continuum 
covariant derivatives by difference operators. Furthermore, as we will argue later, we are naturally 
led to use covariant versions of forward and backward difference operators. The results of these 
difference operators when acting on lattice fields should also result in new lattice fields with well- 
defined gauge transformation properties. Thus hitting a site field with such a (forward) difference 
operator should lead to a field which gauge transforms like a link field etc. It is not hard to show 
that the following definition of a forward difference operator acting on site and link field satisfy 



these properties [ |34| ] 



Dj/vto = U^Mx + ^-Mxp^x+v) (4.6) 

They clearly reduce to continuum derivatives as a — > and ensure that the resulting fields gauge 
transform in the appropriate way. Rather remarkably it is possible to write a covariant backward 
difference which is truly adjoint to the above operator for gauge invariant quantities. Notice that it 
will necessarily act as a divergence operation (for a proof of this statement we refer the reader to 



the paper by Aratyn et al [34]) 



Dpfr to = U (*) U l to -Ul(x-n)fn(x-n) (4-7) 
D^Uv(x) =/ A( vto^ t (^ + v)-^(x-At)/ M v(^- J u) (4-8) 

We have now almost all the ingredients needed to construct our lattice super Yang-Mills theory. The 



remaining key observation goes back to work by Beefier, Joos, Rabin and others [ J35| , |36j] concerning 
discretizing actions written in terms of p-forms and exterior derivatives. They give a topological 
proof using results from homology theory that such actions can be discretized in such a way that 
the resulting lattice theories are completely free of spectrum doubling - that is, the resulting lattice 
theories do not contain any modes having no analog in the continuum theory. However, such 
theories, like staggered fermion models, will necessarily describe degenerate flavors of fermion. 
While this is a problem for QCD it is a bonus in the supersymmetric theories studied here as these 
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degenerate flavors can represent the multiple fermions inherent in extended supersymmetry. While 
these original papers showed how to do this for ordinary derivatives, a generalization to covariant 
derivatives was first given by Aratyn at al.[34] which is the one employed here. The discretization 
prescription is given explicitly by replacing the usual partial derivatives in the continuum theory by 
appropriate difference operators: 

D M -> D+ if it acts like d 
D M -> D if it acts like <i f 

Following this rule we find that the Yang-Mills field strength is discretized as follows 

F^ v {x) = D+U v (x) F^° nt as a - (4.9) 

Notice it is automatically antisymmetric in its indices. Using these ingredients we propose the 
following lattice action for jV = 2 super Yang-Mills 



Sl = ^eTr£(--l)'[0,0]-2^ 2 F 12 (410) 



- ^ 2 B 12 -v^O+?+h.c) (4.11) 

Notice the necessity of introducing both fields and their hermitian conjugates in order to write down 
gauge invariant expressions. The corresponding lattice supersymmetry is a simple generalization 
of the continuum ones [28] 

QUn = V/i 

Q(j) = 
QXn = B n 
QB l2 = [f*i 2 ] (12) 

G? = r] 

Qn = [0,?] (4-12) 

where I have replaced D^ip by D+ip as required by gauge invariance and the prime reflects the fact 
that we must use a shifted commutator as previously discussed. Carrying out the Q-variation and 
then integrating out the field B as in the continuum we obtain the on-shell lattice action 

Sl = fTr£Q[0,0] 2 + 4F 12 (4.13) 
1 

4" 

+ (D+0) t D+0-2 Zl t 2 (D> 2 -D+^ 1 ) (4.15) 
- 2yJfl+!+h.c) (4.16) 
This action is invariant under Q, finite gauge transformations and an additional SO(l, 1) symmetry 

Vn^XWn, n,Xn^jn,Xu (4-17) 



j^l^-xU^Xn^ + vli^V^ (4-14) 
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0^A 2 (4.18) 

This symmetry is the analog of the usual exact U(l) chiral symmetry of staggered quarks. It is 
important to realize however that our coupling of gauge fields to lattice p-form fields does not 
correspond to the usual gauging of staggered quarks. 

It is illuminating to examine the gauge action in more detail. Substituting the definitions of 
Fpy we obtain 

^^ 2 {x)F n {x) (4.19) 

X 

j3Tr£ (21 -U P - U$\ + j3Tr£ (M 12 +M 21 - 21) (4.20) 

X X 

where 

Up = U l {x)U 2 {x+l)Ul{x + 2)ul{x) (4.21) 
resembles the usual Wilson plaquette term and 

M 12 = U l (x)UJ(x)U^(x + l)U 2 (x+l) (4.22) 

is a new zero area Wilson loop term which would vanish if the link variables were restricted to 
unitary matrices. Notice the appearance of the Wilson term depends crucially on the appearance 
of both F^ v and F^ v which lends some support to the use of complex variables in the formulation 
of the theory. Such a theory does not exhibit doubling of course and we might guess that the exact 
supersymmetry will thus prevent doubles from occurring in the fermionic sector also. We turn to 
this question now. 

We can assemble the various p-form components of the Dirac-Kahler fermion back into a 4 
component object in the following way 

fn/2\ 



(4.23) 



\ V2 J 

The twisted fermion action can then be written in the matrix form *F'M*P where 



-AT* 



and the operator 



is just a discrete kinetic term for the fermions. The Yukawas appear on the diagonal as indicated. 
After integration we expect that the result will be a Pfaffian of the operator M (the reality conditions 
in the continuum are important here). To examine potential doubling problems in this operator it 
suffices to examine M in the free limit where its determinant reduces to the usual determinant of a 
two dimensional double free scalar laplacian. 

Pf(M) = det(£) = det(D+D A( ) (4.26) 
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Thus, as advertised, the use of Dirac-Kahler fermions avoids the appearance of spurious zeroes in 
the fermion operator. 

5. Continuum Limit 

The lattice formulation we have given requires a doubling of degrees of freedom - we have 
argued that this is quite natural in any lattice theory and can be associated with the two possible ori- 
entations of the underlying p-cube. We have parametrized this doubling in terms of complex fields. 
However, the target continuum theory that we are hoping to reproduce in the limit of vanishing 
lattice spacing corresponds to setting 

JmX° = all fields X bar scalars (5. 1) 

? = (5-2) 

At this point we must address how this continuum theory can be recovered from our complexified 
theory. First, consider the effective action that results from integrating out the grassmann variables. 
It is clearly gauge invariant if we ignore possible phase problems and replace the Pfaffian by a 
square root of the (gauge invariant) determinant. It also clearly has the right classical continuum 
limit. The only remaining question is whether the Ward identities corresponding to the twisted 
supersymmetry still hold (or more conservatively, hold in the limit of vanishing lattice spacing 
with no additional fine tuning) We conjecture that this may be so and have followed this approach 
so far in our numerical work. 

It is possible to make some progress in understanding why this might indeed be true. First 
parametrize the general GL(N,C) link field U^(x) in terms of a unitary component u ll (x) and a 
positive definite hermitian piece (x) in the following way 

U tl (x)=R li (x)u li (x) (5.3) 

Now consider some Ward identity corresponding to the expectation value of the Q-variation of 
some operator < QO >. It is not hard to show that this object can be computed exactly in the 
semi-classical limit j3 — > °°. The argument follows by recognizing that 

4k<QO>=<Q(OS)>=0 (5.4) 

where Q-exactness of the action has been used. Thus the expectation value is independent of 
coupling P and hence can be computed in the semi-classical limit. This is a standard result of 
topological field theory Now consider the second term in the gauge action 

j3Tr£(M 12 +M 21 -2/) (5.5) 

X 

where 

M12 = Ui{x)ul{x)Ul{x+ l)U 2 (x+l) (5.6) 
and insert the general decomposition of the gauge link given above. The result for M\2 is 

M 12 =R\{x)R\(x+\) (5.7) 
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with a similar result for Mi\. Notice it is independent of the unitary piece u^(x). Thus in the limit 
j3 —y oo each (x) is driven to the identity and the Ward identity can be computed in the space of 
purely unitary gauge links. Thus the Ward identities of the truncated theory should hold at least for 
large /3 since they can be viewed as coming from the complexified theory in the limit of infinite /3 . 
Of course large j3 also corresponds to the limit of vanishing lattice spacing and we see that this line 
of reasoning constitutes an argument that the Ward identities in the truncated lattice theory should 
be realized without fine tuning in the continuum limit. As we shall show in the next section our 
preliminary numerical results are consistent with this. 



6. Simulations 

For our simulations we have taken the gauge links to be unitary matrices taking values in 
U (2). As in the continuum theory the scalars (j> and are taken to be complex conjugates of each 
other 2 . To regulate possible IR divergences associated with the flat directions where [0,0] = we 
have added a scalar mass term m 2 (j)^(j). This will break supersymmetry softly. At the end of the 
calculation we should like to send m — > to recover the correct target theory. This bosonic action 
is real positive semi-definite, gauge invariant and clearly has the correct naive continuum limit. 
To this gauge and scalar action we should add the effective action gotten by integrating over the 
grassmann valued fields. We will represent this as 

det'(M([/,0) t M([/,0)+m 2 ) (6.1) 

where M is the lattice Dirac-Kahler operator introduced earlier. Notice that we have added a gluino 
mass term for the fermions which helps regulate potential zero modes. It is chosen to be identical 
in magnitude to the scalar mass. This ensures the scalar and fermion determinants will cancel for 
zero gauge coupling. In these preliminary simulations we have neglected the issue of the phase of 
the fermion determinant which is reflected in the appearance of just M^M in the effective action. 
Finally, the fourth root reflects the Majorana nature of the continuum fermions. Notice that the 
determinant is also gauge invariant. 

To simulate this model we have used the RHMC algorithm developed by Clark and Kennedy 



[37]. The first step of this algorithm replaces the effective action by an integration over auxiliary 



commuting pseudofermion fields F, F 1 " in the following way 

detl (M*M + m 2 ) = J DFDF^ e ff (M^M+m 2 )~ * f (6 2) 

The key idea of RHMC is to use an optimal (in the minimax sense) rational approximation to this 
inverse fractional power. In practice it is important to use a partial fraction representation of this 
rational approximation 

1 N a 

4^ + 1^73 (6-3) 



2 the antihermitian nature of our basis T a actually ensures that i 
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The numbers a,, /3, for i= 1 . . . N can be computed offline using the remez algorithm 3 . The resulting 
pseudofermion action becomes 

N a 

a F' l F + V F 1 " — '— — —F (6.4) 

It is thus just a sum of standard 2 flavor pseudofermion actions with varying amplitudes and mass 
parameters. In principle this pseudofermion action can now be used in a conventional HMC algo- 



rithm to yield an exact simulation of the original effective action [38]. The final trick needed to 
render this approach feasible is to utilize a multi-mass solver to solve all N sparse linear systems 
simultaneously and with a computational cost determined primarily by the smallest shift /$,-. We use 



a multi-mass version of the usual conjugate gradient algorithm [|39|]. In practice for the preliminary 
simulations shown here we have used N = 12 and an approximation that gives an absolute bound 
on the relative error of 10~ 6 over the entire range of the spectrum (the lower eigenvalue is given 
by m 2 and the upper can be determined from trial simulations using a Lanczos method). Further 



details will be presented elsewhere [ 4C ] 



The numerical results we present here come from simulations where the lattice length takes 
the values L = 2, 4, 8 and the mass parameter m = 0. 1 , 0.05 over a range of coupling j3 = 0.5 — > 4.0. 
Notice that these lattices are roughly equivalent to L = 4, 8, 16 for staggered quark simulations. We 



are currently extending these calculations to larger lattices and smaller mass [4C]. 

We have devoted our initial runs to a check of the simplest possible supersymmetric Ward 
identity - namely the value of the total action. If the latter conesponds to a Q-exact observable 
this should vanish. This fact, together with the quadratic nature of the fermion action allows us to 
compute the bosonic (gauge plus scalar) action exactly and for all values of the coupling constant 
j3 we find 

j8 < S B >= ^RL 2 (6.5) 

where R = 4 is the number of generators of U (2). The following graph shows a plot of ^7 < 5g > 
for a range of couplings j8 and three different lattice sizes using a mass m = 0.1. The bold lines 
shows the analytic prediction based on supersymmetry (for clarity we have added multiples of 
—0.25 to the curves and lines to split up the data from different lattice sizes) While there are 
clearly deviations of order 4-5% at small coupling these appear to disappear at large j8 in line with 
our expectations. 

We have also examined the gauge action and the string tension determined by the smallest 
creutz ratio a = — lnW(l, 1)W(2,2)/(W(1,2)W(2, 1)) and plot these as a function of gauge cou- 
pling in the following plots (m = 0.1 and L = 4). These data are consistent with the existence of a 
confining phase. 

7. Twisting Jf = 4 SYM in D=4 

The JV = 4 super Yang-Mills theory in four dimensions contains four Majorana spinors x ¥ I a 
transforming under the Euclidean Lorentz group SO(4) and an additional SOr(4) R-symmetry 



many thanks to Mike Clark for providing us with a copy of his remez code 
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Figure 1: < S B > 



Figure 2: Gauge action for U (2) model 



group. The Kahler-Dirac twist (see [31], [32]) consists of decomposing those fields into represen- 
tations of the diagonal subgroup 

SO(4)' = diag(S<9(4) x SO(4) R ) (7.1) 

As for two dimensions this twisting procedure is equivalent to regarding the supercharges and 
fermions as matrices: 

¥ = ?]/+ y^iYix + 2iXtivYn7v + (7-2) 

+ ^^va7^7v7a + ^^vipY^YvYiYp (7-3) 

To match the 16 real supercharges of the continuum theory it is natural to take these antisymmetric 
tensor fields real. They can then be regarded as the 16 grassmann components of a single real 
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Figure 3: String tension for U (2) model 



Kahler-Dirac field just as for two dimensions. There will now be a scalar nilpotent supercharge 4 and 
corresponding bosonic superpartners (<jt,A^,B^ v , Wn V x>CuvXp) I n tne en ^ B an d C will turn out to 
be multiplier fields and are integrated out of the final theory leaving two scalars , (which arises 
as a gauge parameter as in two dimensions), one vector A m u and the 4 independent components 
of Wf+vx- These latter four fields will combine with the two scalars to yield the usual 6 scalars of 
JV = 4 super Yang-Mills. 

The form of the <2-variations of these fields under the scalar supersymmetry can be obtained 
in analogy to the situation in two dimensions. 

& = n Qn = [0,?] 

QBp V = [^>,Xnv] Qliiv=B liV 

Q<!> = (7.4) 

In deciding how to write down these transformations one is guided by the identification of physical 
fields. If a field is ultimately to be identified as a propagating field it is necessary to make its Q- 
variation yield another field. Multiplier fields vary into commutators. Notice again that the square 
of the supercharge is an infinitessimal gauge transformation parametrized by the field (j>. 

We hypothesize that, following on from our general arguments, it should be possible to write 
the action as Q-exact function 

S = pQA (7.5) 
where the correct gauge fermion of jY = 4 super Yang-Mills turns out to be 



A 



d 4 xTr 



1 



l[W^p,W vX p} 



the full twisted superalgebra was given by D'Adda et al inQl 
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+ D x W Xflv ) 

+ V^O^+-7]^,^] + -0 ;[1V a[W^ vX , 



1 



V2DfrW vXp] + -C, 



fivkp 



(7.6) 



Several of these terms are in common with the gauge fermion of jY = 2 super Yang-Mills theory 
in two dimensions. The new ones involve the 3 and 4-form fields. Of these the terms involving 
derivatives must be present to generate the correct Dirac-Kahler action for the twisted fermions 
(and will simultaneously generate the appropriate kinetic terms for the W-field). The commutator 
term involving Wn vX coupled to %iiv will generate a quartic potential for the W-field analogous 
to that generated for the scalars and <p . This will allow contact to be made eventually with the 
supersymmetric theory where one expects the scalars and W-field to play similar roles. In the 
same way the commutator term involving and W will also generate the necessary mixed quartic 
couplings between the scalars and the W-field. Carrying out the <2-variation leads to the following 
action 

S = p(S B + S F + S Y ) (7.7) 
where the piece of the action Sb involving the bosonic fields takes the form 



S B 



d A xTv 



B, 



liv \*p.v 



\[W^ p ,W vXp ] +D x W llxv + -B^ 



D^Drf + 1 [$M 2 - ^ [0, WW] [0 , W„ vX ] 



1 



V^D|> W vXp] + -C^ vXp 



and the fermion kinetic terms are given by Sp with 



d 4 xTr 



-*Mv°[/x ¥v]-XnvD x d x 



V2 

^D^Yfi ~ -Jf K nv^p D [p. #vAp] 



(7.8) 



(7.9) 



and Sy contains the Yukawa couplings 

1 



S Y 



d xTr 



1 



1 1 

24! 



^ vX p[^>,K^ vXp ]--Xnv[<l>,Xnv] 



+ V^[0,V^i] + 3|^va[0,^vx] 

1 r 1 V2 

+ X^v[0^ X p,W vX p]-Xnv[¥x,W x ^ v }} 



(7.10) 



Integrating over the multiplier fields B^ v and C flvXp and subsequently utilizing the Bianchi identity 
leads to a new bosonic action of the form 



S B = d 4 xTr 




1 K Xp , W vXp }) + {D X W Xllv f + ^ (Zfy W vXp] f 



[<t>,W^ x }^,W^ vX ] 



(7.11) 
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At this point it is useful to trade the W, 6, and K fields for new variables which will allow contact 
to be made between this theory and one of the conventional twists of .jV = 4 super Yang-Mills. We 
write 



WW 



£jiv?LpT] 



(7.12) 



In terms of these variables and after an additional integration by parts the bosonic action reads 
(notice that terms linear in F^ v cancel) 



Sb 



d 4 xTr 



1 



2 F > 

-r 1,, -r 



1 



(7.13) 



Making the additional rescalings % — > 1% and r] — > 4=7] we find the fermion kinetic term takes the 
form 



Sp = J d 4 xTr 



(7.14) 



where = \s-^ v xpXxp * s the dual field. In these variables the Yukawa's take on the more sym- 
metrical form 



J 4 *Tr 



(7.15) 



The new action can be recognized as the twist of jV = 4 super Yang-Mills due to Marcus [|4l|]. 

Finally we will show how this twisted model may be reinterpreted in terms of the usual formu- 
lation of ,jV = 4 super Yang-Mills theory. First, concentrate on the bosonic action and introduce 
the new fields ((j> = (j>\ + ifo) 

=V[ 1 [1 = 0... 3 



01 

02 



(7.16) 



Then the bosonic action may be trivially rewritten as 



(7.17) 



Notice this is real, positive semidefinite on account of the antihermitian basis for the fields. It is 
also precisely the bosonic sector of the jY = 4 super Yang-Mills action with X' the usual 6 real 
scalars of that theory p2|]. 
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Next let us turn our attention to the fermion kinetic term. From our previous discussion it 
should be clear that the twisted fermion kinetic term is nothing more than the component expansion 
of a Dirac-Kahler action: 



(7.18) 



where *F is the Dirac-Kahler field defined earlier with t] — ► T\/2, 6^ v x — > E^vXpWp m & K ^vXp ~^ 
Naively such an action appears to describe a theory with four Dirac spinor fields - rather 
than the two one would have expected for jV = 4 super Yang-Mills. However it is evident that *P 
obeys a reality condition if its associated Dirac-Kahler field is real. This reduces the action to that 
of two degenerate Dirac fermions. To see this in detail let us adopt a Euclidean chiral basis for the 
/-matrices 

X l\ _ ( -ia\ 
I J Yi ~ \ id J 



(7.19) 



It is straightforward to see that the y matrices (and any products of them) obey a reality condition 

f li =CY»C- 1 (7.20) 

where the matrix C is given by 

With real p-form coefficients this implies a reality condition on *F itself 

W*=C X ¥C 1 (7.22) 

This in turn implies that *F + = C^ T C 1 which makes it clear that the result of integrating over the 
Dirac-Kahler field *F should be the Pfaffian of the Dirac-Kahler operator (let us neglect the Yukawa 
couplings for the moment). This, in turn, will correspond to the product of two Dirac determinants. 

To see this in more detail one can use the reality condition to show that successive columns 
»JA«) f »j/ ^ no t independent but are charge conjugates of each other 



vj/(2) = c( vj/U) ) ip(4) =C W 



(7.23) 



These conditions allow us to rewrite the twisted fermion kinetic term in the conventional form 



1 



(7.24) 



a=l,2 



where the spinors are read off as the first and third columns of the *P matrix in this chiral basis 
/ £-I + 2& + 3 \ ( (Vo-¥o) + i (~¥3 + ¥3)\ 



A, 



2 2 1 ~ l A03 

~~ 2# 02 + 2z'Zoi 
¥o + ¥o) + i (¥3 + ¥i 
\-(Y 2 + ¥2) + i (¥i + ¥i) J 



(¥2-¥ 2 ) + i (¥i-¥i] 



V 



2 + 2 2i '%3 

2^02 ~~ 2i% 0l 



(7.25) 



Here, ^ ± = ^ {x±X*) we the usual self-dual and antiself-dual parts of the field. The Yukawa's 
can also be put in the general form 



- £ XlcfV[x\X a ) 



(7.26) 



a=l,2 
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where P = {/, 75 , 7^ 75 , /I = . . . 4 } and the coefficients Cf are just numbers. This is just the struc- 
ture expected of JV = 4 super Yang-Mills. Notice that both Dirac operators M a including the 
Yukawa interactions possess the symmetry 

(D a )* =CD a C 1 (7.27) 

which shows that their eigenvalues come in complex conjugate pairs and hence the associated 
determinants are positive definite. 

8. Lattice JV = 4 super Yang-Mills 

The lattice action is obtained by discretization of the Q-exact action given in eqn. f7i| Notice 
it is important not to use the same dual variables as Marcus in discretizing the theory as the duality 
transformation is not compatible with the lattice gauge transformations we have described. With 



that proviso the prescription is essentially the same as for the two dimensional theory [31J. We list 
it now for a D-dimensional theory. 

• A continuum p-form field f^.,.^ (x) will be mapped to a corresponding lattice field associ- 
ated with the /^-dimensional hypercube at lattice site x spanned by the (positively directed) 
unit vectors . . . jip} . 

• Such a lattice field is taken to transform under gauge transformations in the following way 

- G{x)f^Xx)G-\x + e^ p ) (8.1) 
where the vector = £y=i Hj. 

• To construct gauge invariant quantities we will need to introduce both f^...^ p and its hermi- 

fi...» p {x) - G(x+e^ p )f^ p (x)G-\x) (8.2) 



tian conjugate /!...» (x). The latter transforms as 



Since for all fields bar the gauge field we will assume f^...^, = La/jSj n p T a where T a are 
antihermitian generators of U (N), this necessitates taking the fields ^ to be complex. 

For a continuum gauge field we introduce lattice link fields U^(x) = e^M = e A ^ T with 
A a ^(x) complex together with its conjugate = e'M*). 

A covariant forward difference operator can be defined which acts on a field f^...^ (x) as 
follows 

D + ^f ll[ ... pLp {x) = U fl (x)f fll .,.n p (x+ii) -f lll ... ilp {x)U }l (x+e lll ... llp ) (8.3) 

This operator acts like a lattice exterior derivative with respect to gauge transformations in 
mapping a p-form lattice field to a (p + l)-form lattice field. 

Similarly we can define an adjoint operator D~ whose action on some field f^...^, is given 
by 

^/ Ml ... M ^^)=/ Ml ...^(x)^(x+e Ml ... Mp -M)-^(^-M)//xi...M P (^-M) (8-4) 
It thus acts like the adjoint of the exterior derivative. 
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• As for two dimensions all instances of dp in the continuum action will be replaced by 
D+ if the derivative acts like d (curl-like operation) and D~ if the derivative acts like d^ 
(divergence-like operation). 

This prescription preserves the twisted scalar supersymmetry and gauge invariance and targets a 
complexified version of the jV = 4 super Yang-Mills model. Again, to obtain jV = 4 SYM we 
conjecture that the Ward identities corresponding to the scalar supersymmetry will be satisfied 
without fine tuning in the continuum limit /3 — > °° for the truncated lattice model in which the fields 
are taken real. 

9. Conclusions 

In this paper we have presented a new way to put supersymmetric theories on the lattice. The 
approach only works for theories with a number of supercharges which is a integer multiple of 
2 D if D is the dimension of (Euclidean) spacetime. This includes the important cases of jV = 2 
and jY = 4 super Yang-Mills theory in two and four dimensions respectively. The basic idea 
is to start from a twisted formulation of the supersymmetric theory in which the fermions are 
represented as grassmann valued antisymmetric tensor fields. The original supersymmetry algebra 
is replaced by a twisted algebra which contains a scalar nilpotent supercharge Q. Furthermore the 
action of the theory can then be written as the ^-variation of some function. It is straightforward 
to generalize the action of the scalar supersymmetry to the lattice while maintaining its nilpotent 
property. Invariance of the lattice theory under the scalar supersymmetry can then be obtained 
easily. We have described the details of our discretization prescription in detail. It requires the 
introduction of particular covariant forward and backward difference operators and non standard 
gauge transformation properties of the fields. The final theory is gauge invariant, free of spectrum 
doubling and possesses an exact supersymmetry. The price one pays for this is that the lattice fields 
are in general complexified. We have argued that the Q-exactness of the action ensures nevertheless 
that the theory truncated to the real line will satisfy the Ward identities corresponding to the scalar 
supersymmetry without fine tuning in the continuum limit. We have also presented preliminary 
numerical results stemming from a dynamical fermion simulation of the D = 21] (T) model using 
the new RHMC algorithm. 

There are several directions for further work. The most obvious is the need for both perturba- 
tive and numerical checks of both the exact and broken Ward identities corresponding to the twisted 
supersymmetry. Exactly how much residual fine tuning would be required in the four dimensional 
theory is currently unclear and could be revealed by such an analysis. It should also be possible 
to use these ideas to construct lattice regularizations for the jY = 4,8 theories in two dimensions 
and the .jV = 8 theory in three dimensions. This would allow contact to be made with predictions 
from string theory and supergravity theories in lattice models which could be simulated relatively 
easily and with high precision. The geometrical nature of the twisted theory lends one to speculate 
that it should be possible to formulate these models on general simplicial lattices. If this could be 
realized, then the resulting models could be coupled to gravity in the usual way by carrying out a 
sum over random lattices resulting in a lattice regulated supergravity theory. It would also be very 
interesting to understand the relationship, if any, between twisting and the orbifold constructions 
of Kaplan at al. 
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It appears that the twisted formulation of a theory with extended supersymmetry affords a more 
useful starting point for constructing lattice models over conventional approaches based on spinors. 
However, the precise details of the discretization procedure are not determined uniquely. While this 
work was in preparation a new paper by D'Adda et al [43] was posted which claims to maintain 
all the supercharges of the two dimensional twisted theory by adopting a different mapping of the 
continuum twisted theory to the lattice. This appears to be a very intriguing development. 
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